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C^") , Abstract. The Laser Interferometric Space Antenna (LISA) will observe supermas- 

sive black hole binary mergers with amplitude signal-to-noise ratio of several thou- 
sands. We investigate the extent to which such observations afford high-precision tests 



f**^ , of Einstein's gravity. We show that LISA provides a unique opportunity to probe the 

^? ' non-linear structure of post-Newtonian theory both in the context of general relativity 

£jl and its alternatives. 
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In a black hole binary, as the two holes orbit about their centre-of-mass, the energy 
and angular momentum from the system is dissipated into gravitational radiation. The 
radiation back-reaction force causes the two bodies to gradually spiral in towards each 
other, resulting in a strong burst of radiation just before they merge to form a single 
black hole. The radiation emitted at the end of the binary evolution is the primary 
target of both the ground- and space-based interferometric gravitational wave (GW) 
antennas pQ. It has become necessary to gain an accurate understanding of the late-time 
evolution of binaries in order to help data analysts in detecting the signal and measuring 
the parameters by fitting the observed signal with that expected from general relativity. 

In General Theory of Relativity there is no exact solution to the two-body problem. 
In the absence of an exact solution (analytical or numerical) theorists have resorted to 
an approximate solution to the problem using post-Newtonian (PN) theory. Briefly, the 
programme of PN theory is the following. Let us consider a binary consisting of two 
non-spinning black holes of masses mi and m 2 (total mass M = mi + m 2 , symmetric 
mass ratio rj = mim2/M 2 ) separated by an orbital (Schwarzschild coordinate) distance 
r(t). Post-Newtonian theory expresses the relevant physical quantities as a series in 
the 'small' velocity parameter e = v/c = ^jGM/rc 2 . For a system consisting of non- 
spinning black holes the only relevant quantities are the (specific) binding energy E 
and the GW flux JF, which are obtained as perturbative expansions in v. (Following 
the standard convention, in units where G = c = 1, v n corresponds to a term of order 
7^-PN.) Currently, these expansions are known to order v 7 or 3.5PN J2J0]: 

l _ 3 _ qo _ 7 _ 

E = -- V v 2 Y,E k v 2k , ^=^V°^^ fc , (1) 

fc=0 fc=0 

where Ek and T^ are the PN expansion coefficients that are all functions of the mass 
ratio rj. 

In the adiabatic approximation one then uses the energy balance equation, 
—dE/dt = J 7 , to compute the evolution of the orbital phase tp(t) using the following 
coupled ordinary differential equations: 

dip v 3 dv dE/dt —T 

~dl =UJ = Jd ) dt = dE/dv = E'{v) ' ^ 

where E'(y) = dE/dv. The phase $(£) of the emitted radiation at dominant order is 
simply twice the orbital phase: $(£) = 2ip(t). 

The phasing formula obtained by solving the above differential equations includes 
different PN terms arising from nonlinear multipole interactions as the wave propagates 
from the source's near- zone to the far-zone [I]. The 1.5PN and 2.5PN term arise solely 
due to the interaction of the Arnowitt-Deser-Misner (ADM) mass of the source and the 
quadrupole moment. It is physically due to the scattering of quadrupolar waves off the 
Schwarzschild curvature generated by the source and is referred to as the gravitational 
wave tail. The 3PN term includes, in addition to the terms at the retarded time, more 
interestingly the cubic nonlinear interactions due to the scattering of the wave tails by 
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the ADM mass energy of the spacetime. The observational tests of these PN terms in 
effect test the nonlinear structure of Einstein's gravity. 

In the restricted PN approximation (so called because one keeps only PN corrections 
to the phase of the radiation but neglects amplitude), the response of an interferometric 
antenna to the incident radiation from a source at a luminosity distance Dl is 

h{t) = -^ [7iMF(t)] 2/3 cos $(£), (3) 

where A4. = rf'^M is the so-called chirp mass of the system, F(t) = -^—^ is the 
instantaneous frequency of the radiation, < C < 1 is a dimensionless geometric factor 
that depends on the relative orientation of the binary and the detector and its average 
over all orientations is C = 2/5. The importance of including PN corrections in the 
phase, and neglecting them in the amplitude, was realized quite early on |3] and has led 
to a lot of simplification of the data analysis problem. For the tests proposed in this 
Letter it may eventually be necessary to incorporate these amplitude corrections [S]. 
For the sake of simplicity, however, we have refrained from doing so in this work. 

For our purposes it will be useful to work with the Fourier transform h(f) = 
J_ h(t) exp(2iiift)dt of the above signal. Using the stationary phase approximation it 
has been shown that [TJ Ej 



h(f)=Ar 7 / 6 ex V ^(f) +l ^\, (4) 

with the Fourier amplitude A and phase ^{f) given by 

A =— A /^A4 5/6 

*(/) = 2nft c + $c + Yl M {k ~ 5)/3 - ( 5 ) 

Here tc and $<7 are the fiducial epoch of merger and the phase of the signal at that 
epoch. The coefficients in the PN expansion of the Fourier phase are given by [Hll^llUj: 

^ = -^(nMf-^a k , (6) 
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Unlike the other parameters ipk, k 7^ 5,6, which depend only on the masses of the system 
and are thus constants, the parameters ^5 an d V'e are n °t constants and have some 
log /-dependence. We treat the log-terms as constants with the justification that the log- 
dependence on the frequency is weak in the relevant bandwidth. Another possible choice 
is where the coefficients of the log-terms are treated as additional signal parameters. This 
choice, indeed, increases the dimensionality of the parameter space making the Fisher 
matrix highly ill-conditioned and we do not discuss it here. 

Post-Newtonian theory has been highly successful in explaining the decay of the 
orbital period in binary pulsars and in confirming the emission of gravitational radiation 
by these relativistic systems (cf. ^T] and references therein). Nevertheless, the binary 
pulsar radio observations do not test PN theory to a high order. This is because the 
typical velocity in the most relativistic of binary pulsars is v ~ 3 x 10 -3 , which is not 
large enough for higher order terms to be important. 

GW observation of the coalescence of binary black holes (BBH) will provide a unique 
opportunity to test the PN theory to very high orders. This is because the velocities 
in the system, close to the merger, could be as high as v ~ 0.2-0.4, making the highest 
order known PN term 10 12 -10 14 times more important for GW observations than for 
radio binary pulsars and several tests of general relativity have already been proposed 
by various authors ^21 CHI EI- We shall show in this Letter that the brightest events 
that can be expected in the space-based Laser Interferometer Space Antenna (LISA) 
will test all the PN terms computed so far. 

Our proposal to test the PN theory is the following: let us suppose we have a GW 
event with a high signal-to-noise ratio (SNR), say more than 1,000. Once an event is 
identified we suggest to fit the data to a signal wherein each term in the PN expansion 
is treated as an independent parameter. More precisely, instead of fitting the detector 
output with a signal that depends on only the two mass parameters, we could fit it with 
the same signal but by treating all the t/Vs as independent. For example, if we want to 
test the PN theory to order v A then we should use a four-dimensional grid of templates 
consisting of {ipo, ip2, ^3, ^4} rather than the two-dimensional one that is used in the 
detection problem. This higher- dimensional fitting of the data with our model would 
measure each of the PN coefficients independently of the others. In Einstein's theory, 
for the case of nonspinning binaries, each of the f/Vs has a specific relationship to the 
masses, ipk — ^A;( m i5 m 2) 5 whereas in a different theory of gravity (for example, a theory 
in which the graviton has non-zero mass) the relationship will be different and might 
involve new parameters. Thus, the measured ^'s could be interpreted, in principle, in 
the context of different theories of gravitation. 

In the case of general relativity we know that the f/Vs are given in terms of the 
masses by Eqs. © and Q. If general relativity (or, more precisely, the PN theory that 
approximates general relativity) correctly describes the dynamics of the system then the 
parameters must be consistent with each other within their respective error bars. One 
way to check the consistency would be to invert the relationships between the ^'s and 
the masses to obtain m 2 = m^mi, ^), and plot m 2 as a function of m 1 for various 
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4>k s i an d see if they all intersect at a common point. If they do, then the theory is 
correct to within the measurement errors, if not, the theory is in trouble. In addition to 
the PN theory we could also test other approximants, such as the P-approximant fH] 
or the effective one-body approximation jTH] that have been proposed as alternatives to 
the orbital dynamics of binary inspirals as also numerical relativity predictions. 

Although these tests are in principle possible, an important question is whether 
the various PN coefficients can really be measured accurately enough for the test to be 
meaningful. We already know that a simpler test proposed in Ref. [I2j, in the context 
of ground based detectors, requires events with SNRs in excess of 25. The generalized 
tests proposed in this Letter would require much stronger signals, SNRs of 100 to test 
lower-order terms and in excess of 1,000 to test all terms currently known. Initial ground- 
based interferometers are unlikely to observe events with such large SNRs. While some 
of these tests might be possible with advanced detectors, a supermassive BBH merger 
in LISA is our best bet. 

To test an approximation it should be possible to measure the various PN 
coefficients with a good accuracy. We shall require that the relative error in the 
measurement of a parameter be less than 100%, i.e. Aipk/ipk < 1, where Aipk is the 
error in the estimation of the parameter -0*., in order that its presence is tested with 
confidence. 

Next, we will work out the accuracy with which we can measure the various 
parameters in the standard PN theory. Given a signal h(f; 9k) that depends on the 
parameters 9k the information matrix g a b is given by 

9km = (hk, h m ) , (8) 

where h m = dh(f; 9k)/d9 m , and for any two functions a and b their scalar product (a, b) 
is defined as: 

(a, b) = 2 r \a(f)b*(f) + aVW)} T^L, (9) 



s h (f) 

with a(/) being the Fourier transform of a(t). The quantity Sh{f) is the one-sided 
noise power spectral density (PSD) of the detector with units of Hz" 1 ; it is the only 
characteristic of the detector that enters the calculation of the errors in the estimation 
of parameters. 

The variance-covariance matrix C a b is nothing but the inverse of the information 
matrix, Cab — 9ab- The diagonal components Ckk are the variances in the parameters 
9k- Ckm I VCkkC mm , k j^ m, are the correlation coefficients between parameters 9k and 
9 m , taking on values in the range [—1, +1]. 

For the binary system in question defining a new parameter vector = 
{t c , $c, ipo, 1P2, ■ ■ ■ , ^7} one finds 

got = 4A 2 / /1S ° f ab r 7/3 



fa 



s h (fV 



_ d9(f;6)d9(f;6) 
U = ~W a d9b~- (10) 
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Note that we have now introduced specific values for the limits in the integration: f s is 
simply the lower frequency cutoff chosen such that the loss in the SNR due to this choice 
is negligible compared to the choice f s = and /i so is the frequency of the radiation 
at the last stable orbit of the system which we assume to be the value given in the 
test mass approximation, namely /i so = l/(6 3//2 7rM). Using the Fourier phase given in 
Eq. (0) it is straightforward to compute f a t : 
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where x = f 1 ^ 3 . We see that the information matrix will involve moments of the noise 
spectrum of the form / °° f~^ 3 s jj\ , where j runs from 1 to 17. The elements of the 
information matrix, therefore, take on values in a very large range leading to a highly 
ill-conditioned matrix. Extreme caution should be exercised in computing the moments, 
else it is easy to end up with values in the covariance matrix that are negative, and even 
imaginary, while we know that the covariance matrix should be real. We have verified 
our results using three different methods for numerical integration and found that they 
differed from one other by no more than 10 %. 

We assumed that LISA consists of only one interferometer with sensitivity as in 
Ref. ^7J E3 [T3j and the binary consists of two black holes of equal masses in quasi- 
circular orbit and observed for one year before merger. Fig. Q plots the error in the 
various parameters ipk as a function of the total mass for optimally oriented supermassive 
BBH mergers at a distance of Dl = 3 Gpc. We shall be interested in binaries with their 



10 



1U 

10 2 

10 1 

*10° 


A V( /V 

N A\|/,/\|/, 

^ cv --Avi< 4 /y 4 


< 10" 1 


"•••••,., _.~~- ,--"'.••• 


10" 2 


"^ 


10" 3 


^ _^^ 


-,„-4 





10° 



10 10 10 

Total Mass of the Binary (M ) 



10' 



Figure 1. Errors in the measurement of the PN coefficients. The sources are optimally 
oriented supermassive BBH located at a luminosity distance of Dl = 3 Gpc. 
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total mass in the range 2xlO 4 M — 2xl0 7 M Q . At 3 Gpcsuch sources would produce SNR 
in the range 1,000-6,000. From Fig. Q wherein the log-terms are treated as constants, 
it is clear that in the mass range 10 5 — 10 6 , fractional errors associated with most of the 
parameters are less than 1. Thus, LISA will provide a unique opportunity to test the 
PN and related approximations to a high degree of accuracy. 

Fig. ^illustrates how one might test for the consistency of the individual masses for 
the specific system of the merger of black holes each of mass 1O 51 M . We have chosen 
a system that gives the lowest errors in different parameters (cf. Fig Q) and assumed 
that the source is at a distance of 1 Gpc. For each parameter ipj, we have plotted the 
region enclosed by the boundaries ipk + Aipk and ipk — Aipk, where Aipk is the one-sigma 
error in the estimation of ipkt- The higher order ^.'s will have to enclose the region 
determined by, say ipo and ip2- This will be a stringent test for the various parameters 
and will be a powerful test if LISA sees a merger event with a high SNR of ~ 10 4 . For 
binaries that merge within 1 Gpc the test would confirm the values of the known PN 
coefficients to within a fractional accuracy of ~ 1. 

It would be interesting to ask whether the proposed tests can, in principle, 
distinguish general relativity (GR) from, say a theory that also includes a massive 
graviton [T^]. In this theory the 1PN parameter ip 2 is different from that in GR. 
Remarkably, we find that an year's worth of observation of BBH mergers in the mass 
range 2 x 10 4 — 2 x 10 7 M Q should be sufficient to discriminate GR from a massive graviton 
theory provided the Compton wavelength of graviton \ g < 5.5 x 10 14 — 3.8 x 10 15 kms. 
These limits make the simplifying assumption of neglecting the as yet uncomputed 

\ Since tpe has only log-dependence on the frequency, we have treated it to be a constant and have 
chosen / = /i so in -06 in our calculations. Other choices of the frequency would alter the shape of the 
contour corresponding to -06 in the toi — m^ plane, and our choice corresponds to the one for which i/jq 
will be the largest as the signal sweeps the detector band and leads to smallest errors. Finally, there 
is no test corresponding to ip$ in Fig ^ an d El as a constant i/>5 gets absorbed into a redefinition of the 
coalescence phase $c- 
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Figure 2. The mi-m,2 plane plot for different f/Vs for the merger of a BBH of 
2 x 10 51 Mq, at a luminosity distance of lGpc. 
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higher PN order corrections to GW phase in the massive graviton case. (See also ^H] 
for a discussion regarding the extent to which GW observations can critically distinguish 
between different theories of gravitation in comparison to the binary pulsar tests.) 

In the present work we have dealt with only non-spinning binaries. The spin 
parameters, (3 from spin-orbit coupling at 1.5PN and a from spin-spin coupling at 2PN, 
are assumed to be less significant for these equal mass systems. For the unequal mass 
case the spin effects are expected to be more important. Orbital eccentricity, which 
might introduce systematic effects in these tests, has not been included as we have 
restricted our analysis to binaries in quasi-circular orbits. Massive graviton theories can 
be tested since they lead to a phasing formula that is structurally similar to general 
relativity but with terms modified due to the propagation delay. Alternative theories 
of gravity, such as the Brans-Dicke theory, where the PN structure of the phasing is 
different due to the presence of dipolar radiation, may also be tested by a straightforward 
extension of the above proposal. These and other issues will be investigated in future. 
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